Increasing thermoelectric performance using coherent transport 
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We show that coherent electron transport through zero-dimensional systems can be used to tailor 
the shape of the system's transmission function. This quantum-engineering approach can be used 
to enhance the performance of quantum dots or molecules in thermal-to-electric power conversion. 
Specifically, we show that electron interference in a two-level system can substantially improve 
the maximum thermoelectric power and the efficiency at maximum power by suppressing parasitic 
charge flow near the Fermi energy, and by reducing electronic heat conduction. We discuss possible 
realizations of this approach in molecular junctions or quantum dots. 
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Thermoelectric devices are currently of high interest 
both for solid-state cooling and for increasing energy ef- 
ficiency in converting heat into electric power. The effi- 
ciency of thermoelectric materials to convert a tempera- 
ture gradient into electrical work is characterized by the 
figure of merit ^ 



ZT = 



S^aT 



(1) 



with the temperature T, the Seebeck coefficient 5, the 
electric conductance a, as well as the electronic and 
phononic thermal conductance tZei and hZph- Other im- 
portant performance parameters are the power output 
P = —{jiL — fJ^RjI^ the maximum power Pmax at the 
optimal bias and level configuration, and the efficiency 
T] = P/Jq. Here, /i^ and jj.r are the chemical poten- 
tials of left and right bath, respectively, / is the particle 
current, and Jq is the heat flux out of the warm bath. 

Nanostructured materials are attractive candidates for 
efficient thermoelectrics, because they offer the oppor- 
tunity to optimize S by using the energy-selectivity 
of charge carrier transport in low-dimensional systems^ 
combined with the suppression of hiph by interface 
scattering^ In such systems, /, S^ a, Jq and Kqi can be 
directly evaluated from the electronic transmission func- 
tion ^{E).^ 

In particular, zero-dimensional (OD) systems, such as 
quantum dots or molecules, that are weakly coupled 
to electron reservoirs, can be designed as ideal energy 
filters with an energy-dependent transmission function 
Ti{E) oc S{E — £^i), where Ei is the position of the sin- 
gle level that contributes to transport. In this limit, and 
for hZph = 0, ZT diverges^ and the efficiency of ther- 
moelectric power conversion approaches Carnot efficiency 
ivc)-^'^ Additionally, such ideal quantum dots reach the 
ideal Curzon-Ahlborn efficiency r]cA of about 7^c/2 when 
El is tuned to maximum power production Pmax-^ How- 
ever, the limit ^{E) oc 5{E — Ei) in OD systems is not 
interesting for applications because the power output be- 
comes exceedingly small, and even a small hZph ^ leads 
to a low value of ZT. To increase the current and Pmax in 
OD systems one needs to broaden ^{E). This drastically 
reduces the efficiency at maximum power r]maxP^ because 



the Lorentzian-shaped S(^) has a long low-energy tail 
that leads to a counter flow of cold charge carriers, con- 
tributing with opposite sign to S.^ 

Here, we show that coherent transport in OD systems 
can be used to tailor Ti{E) such that counter flow of cold 
charge carriers is effectively suppressed. In addition to 
its fundamental interest, this approach to quantum engi- 
neering is shown to substantially increase Pmax^ VmaxP 
and ZT compared to the ideal OD systems addressed 
above. We discuss how this effect can be implemented in 
semiconductor quantum dots or in molecular junctions. 

We consider a two-level system with both energy lev- 
els Ei^E2 situated on one side of fiL and fiR^ (Fig.[Ha)), 
such that the charging energy is of minor importance^ 
and a spin degeneracy will not affect transport signif- 
icantly. This setup is similar to the double-dot case^^ 
addressed very recently in Ref. luJ, where equal coupling 
strengths of the two levels were assumed. Related find- 
ings were presented in Ref. [l2| where the thermoelec- 
tric properties of double quantum dots with couplings 
ranging from serial to parallel configurations were inves- 
tigated, assuming equal energies of the two dots. 

The coupling between lead and dot is parametrized as 
in Ref. 9, assuming equal coupling strengths to left and 
right leads. The two levels couple with different parity to 
the leads and their coupling strengths differ by a factor 
a^, Fi^i = Fi^i = F, Tl2 = ^R2 = ci^^- This difference 
in parity will turn out to be essential for the increased 
thermoelectric performance. For such couplings it can 
be shown that the Breit-Wigner formula provides results 
identical to the exact method of nonequilibrium Green's 
functions (NEGF)!^ 



E(^) 



E- Ei^iT E - E2^ io?T 



(2) 



where the two levels contribute with different sign due 
to the opposite parity. Vanishing transport at the Fermi 
energy Ep = corresponds to E2 = o?Ei. A device 
operated close to such a level configuration would have 
the desired property that S(0) ~ independently of F. 

Fig. [2] shows that a large power factor S'^cj is achieved 
for level positions over an energy range as large as several 
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FIG. 1. (a) Level positions in the proposed two- level system. 
(b)-(c) Possible implementation using zinc porphine, sketch- 
ing the wave functions of the two almost degenerate levels 
close to the Fermi energy, see also Ref. 14. The light /dark 
shaded regions indicate a positive/negative amplitude of the 
wave function ^, respectively. By replacing the hydrogen 
atoms at the positions L and R with, for example, sulfur the 
molecule can be contacted to leads so that the parity of the 
wave functions with respect to the leads differ for the two 
levels. 



ksT, {ksT ^ 0.26 eV at T = 300 K). The transmission 
function at maximum power factor, which is marked by 
a cross in Fig. [21 is displayed in Fig. [3fa). The strong 
asymmetry in S(£^), with a sharp step facing Ep^ is ideal 
for high power production. ^^ The decrease of ^{E) for 
large energies prevents the transmission of high-energy 
electrons, which reduces tvel and results in an increased 
efficiency of the device. 
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FIG. 2. The power factor for coherent transport through a 
quantum system with two levels and for T = 300 K, F = 
SksTn, and a — 0.6. Level positions are defined relative to 
Ef' The cross marks the maximum (S^a — 316 fW at Ei — 
5.80kBT, E2 = 1.47/cbT), and the transmission function is 
calculated at this level configuration in Fig. [3ja). 

Fig. [3l^b) displays Pmax(r) for different values of the 
asymmetry parameter a, showing peak power values 
more than two times larger than the Pmax achievable in 
a single- level OD system corresponding to a = 0. In the 
coherent case maximum power occurs for larger F, since 
the broadening of the levels is not such a severe problem. 
Fig.[3](c) shows the efficiency at maximum power r]maxP^ 
which can be increased by over an order of magnitude, 
compared to the single level case, for large F. In the limit 
F ^ one observes that r]maxP -^ Vca ^ 0.br]c-^^^ When 
F and a are small Pmax and r]maxP are not substantially 
increased, since the effect of the second level is small 
when ksT ^ a^T. As F is increased, the effect of the 
second level can be observed as an increase in Pmax and 
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FIG. 3. (a) Transmission function at the level configuration 
corresponding to maximum power, see Fig. [2] (b) Maxi- 
mum power production and (c) efficiency at maximum power 
rjmaxP in units of the Carnot efficiency 77c, as a function of 
F. (d) shows rj^axP (defined as Pmax at the optimal choice 
of F, see text), with the dashed line marking the one- level 
result. In (b)-(d) the temperatures of the two leads are given 
by Tl = 330 K and Tr = 300 K. 



VmaxP^ which can result in a local maximum. Fig. ^d) 
displays Tj^axP^ ^^^ efficiency at maximum power where 
the power production is optimized with respect to F as 
well as with respect to bias and level positions. We re- 
strict ourselves to F < lOksTR where effects of the local 
maximum, present for small a, is not observed. For small 
a the presence of the second level is negligible, and the 
efficiency of a single level is approached. High efficiency 
is achieved close to E2 = o?Ei. For a ^ 1 the level 
configuration that yields Pmax does not coincide with 
E2 = o?Ei where transport is blocked (Eq. [2]). Thus 
VmaxP decreases for large a and the maximum is found 
around a = 0.7. 

To facilitate the comparison with other work we 
present, in Fig. HJ^a), the resulting figure of merit ZTei, 
defined by hZph = in Eq. ([1]). We see that even for 
the relatively large F resulting in maximum power, val- 
ues as high as ZT^i = 7 can be reached. In part, this 
increase is due to a decrease in Kqi. It is worth noting 
that, for the realistic case of finite Hiph, maximum ZT 
is expected near the conditions for Pmax (indicated by 
crosses in Figs. [2] and HJ^a)), because a high power fac- 
tor is needed in Eq. ^ to provide robustness against 
parasitic heat flow due to phonons. 

Coherent transport through molecules has been dis- 
cussed previously in Refs. [l6| andllTl, where the interfer- 
ence between levels positioned on different sides of Ep 
causes a dip in the transmission function at Ep. The 
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FIG. 4. (a) Thermoelectric figure of merit ZT^i for cofierent 
transport as a function of level configuration for T — 300 K. 
The parameters V — SksT, and a = 0.6 match the conditions 
for maximum power in Fig. [3{b). (b) The thermopower S 
calculated for [/ = and t/ = 10 /csTmean, F = /csTmean, a = 
0.7, and AT = 0.1 /csTmean. The second level is positioned 

at E2 — kBTmean- 



corresponding reduction of charge flow near Ep can be 
used to increase the thermoelectric figure of merit, for 
example in molecules^^ or quantum dots.^^ However, the 
increased values of ZT^i obtained in these papers are a 
result of the decrease in Kqi. This results in low power 
output and in a high sensitivity to phonon^ because 
ZT = ZTei/{l + i^ph/f^ei)- Increased thermopower due 
to interference effects has also been discussed in the con- 
text of side-group induced Fano resonances. ^^ While the 
thermoelectric performance was significantly improved, 
this approach is challenging because it requires very good 
control (~ 10°) of the tilting angle of the side group. 
Recently, efforts were made to determine which types 
of molecules exhibit Fano resonances, ^^ and investigate 
their thermoelectric performance.^^ 

We now discuss possible experimental implementations 
of coherent enhancement of thermoelectric efficiency: 
First, highly controllable CD systems can be realized in 
the form of quantum dots defined in two-dimensional 
electron gases^^ or in nanowires^i^^ By electrostatically 
gating the quantum dot, its energy levels can be shifted 
so that mainly two states contribute to the transport. 
The relative position of the two levels with different par- 
ity can be controlled by applying a magnetic field, ^^ in 
particular in quantum dot materials with a large g-factor 
such as InSb.^^ The use of two parallel coupled quantum 
dots allows for better control of the two levels, but co- 
herence is decreased due to the spatial separation of the 
two dots. 

The exquisite experimental control over semiconduc- 
tor quantum dots thus makes them an excellent testbed 
for the fundamental study of the effects predicted here. 
However, the use of quantum dots in applications, for 
example in the form of nanocrystals embedded into bulk 
materials, is limited because it is very challenging to pro- 
duce such dots with the high uniformity required to have 
equal energy levels in all dots. 

An alternative system with inherently much better re- 
producibility are molecular junctions that fulfill the fol- 



lowing requirements: 

i) The molecule must have two almost degenerate, con- 
ductive states with different parity. Such states can be 
found in molecules that are symmetric with respect to 
the leads, where the symmetry implies that the states 
must be symmetric and anti-symmetric, 
ii) The two states must be positioned within a few F of 
Ef (Fig. E]). The active states could be either LUMO- 
states which should be placed slightly above Ep, or 
HOMO-states which should be placed slightly below. 

A possible candidate is zinc porphine (Fig. [I{b,c)). 
This molecule has two almost degenerate HOMO-levels 
with the desired parity^ The relative position of the two 
levels can be changed by adding e.g. phenyl groups^^' 
or replacing all hydrogen atoms with e.g. chlorine or 
boronJ^ To contact the molecule by gold electrodes one 
could replace two hydrogen atoms by sulfur (Fig.[T]^b,c)). 
The leads can affect the molecular orbitals via charge 
transfer^^ and symmetry breaking. To include such ef- 
fects, and calculate the transmission function of the sys- 
tem, one can use the method described in Ref. 29. A 
rough estimate of the coupling between lead and orbital 
i of the junction molecule is given by F^ ^ 27r|^ipC^A^0 5 
where ^^ is the value of the wave function at the sulfur, 
C = 2 eV is the matrix element between the gold and 
the sulfur atoms, and Nq = 0.07/eV is the gold density 
of states per atom per eV at EpS^ Assuming that the 
wave function is evenly distributed among the 20 carbon, 
4 nitrogen, 1 zinc, and 2 sulfur atoms of zinc porphine 
gives ^i = (27)-i/2 and F^ ^ 0.065 eV ^ 2.5 ksT at 
T = 300 K. The actual values for F^ will vary between 
the different orbitals, but this estimate shows that the 
parameters are in a regime where coherent transport can 
be significantly more efficient than transport through a 
single level, (Fig. [3jb) and (c)). 

We now turn to the effect of phonons on thermoelec- 
tric performance. Values around hCph = 50 pW/K have 
been experimentally determined for molecular wires ^ 
and Kph = 10... 100 pW/K has been estimated for molec- 
ular junctions. ^^ For the parameters of Fig.HJa) the lat- 
ter results in impressive ZT = 0.7... 3. These values are 
quite preliminary as the exact value of hZph depends on 
the junction molecule. However, it should be noted that 
experiments and theoretical calculations of the IR spec- 
tra of zinc porphine show that the majority of the vi- 
brational transitions lie significantly above ksT at room 
temperature,^^ suggesting a low hZph- Also any destruc- 
tive effect of phonons on the coherence should not be a 
severe problem as the suggested operation point of the 
device is in the regime F > /c^T, i.e. the coherence time 
of the electron on the junction molecule is longer than 
the tunneling time, resulting in coherent transport. A 
more detailed study of the couplings between the low- 
energy vibrations and the electronic degrees of freedom 
is required. 

So far we have not discussed the effects of electron- 
electron interaction. In the following we show that a 
finite charging energy U does not affect the results qual- 



itatively, and may in fact result in an increased S. To in- 
clude a finite U we use the second order von Neumann ap- 
proach (2vN),^^ an equation of motion technique where 
co-tunneling as well as the coherence and charging energy 
between the levels are included. This approach calculates 
the current without directly calculating the transmission 
function S(£^), and S is calculated as the open-circuit 
voltage at a small but finite temperature gradient AT 
applied between the leads. 

The 2vN-approach neglects certain tunneling events 
of higher order which results in less accurate results for 
lower temperatures. In Fig. HJ^b) the thermopower is cal- 
culated for the temperatures ksTL = /^B^mean + AT/2 
and ksTji = /c^Tmean — AT/2. To investigate the im- 
portance of a finite /7, E2 has been placed slightly above 
Ef at E2 = kBTmean- The thermopower is then calcu- 
lated as a function of Ei . To demonstrate that the 2vN 
method gives accurate results, and that the effect of the 
finite AT is negligible, the U = results are compared to 
results obtained using NEGF. When Ei ^ the trans- 
port is dominated by this level, which results in a low 
thermopower in the case of U = 0. For a finite U part 
of the spectral density of Ei is shifted to Ei -\-U due to 



the small but finite occupation of level E2. This results 
in transport at increased energies, and consequently a 
higher thermopower. When Ei is far away from Ef the 
difference between U = and U = 10 /c^Tmean decreases 
as the effect of the charging energy vanishes when the 
levels are empty. 

In conclusion, we have shown how coherent transport 
can be used to tailor the shape of the system's transmis- 
sion function, to combine a large thermoelectric power 
production with a high efficiency. Compared to transport 
through one- level systems, the maximum power produc- 
tion is more than doubled. At the same time efficiencies 
at maximum power up to r]maxP = 0.33 are achieved 
for f<iph = 0, comparable to those theoretically achiev- 
able in ideal (ballistic) one-dimensional devices where 
VmaxP = 0.36 is expected.^ We also proposed a specific 
molecular junction for the implementation of this novel 
approach to thermoelectric engineering, and expect a ZT 
of the order of unity or above. 
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